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C \ Abstract 
I—* 

■ The spatial correlations in transmitter node locations introduced by common multiple access proto- 

(N ' 

cols makes the analysis of interference, outage, and other related metrics in a wireless network extremely 
difficult. Most works therefore assume that nodes are distributed either as a Poisson point process (PPP) 
or a grid, and utilize the independence properties of the PPP (or the regular structure of the grid) to 
analyze interference, outage and other metrics. But, the independence of node locations makes the PPP 
a dubious model for nontrivial MACs which intentionally introduce correlations, e.g. spatial separation, 
while the grid is too idealized to model real networks. In this paper, we introduce a new technique based 

^ , 

CN| ■ on the factorial moment expansion of functionals of point processes to analyze functions of interference, 

00 

£T} • in particular outage probability. We provide a Taylor-series type expansion of functions of interference, 

wherein increasing the number of terms in the series provides a better approximation at the cost of 
increased complexity of computation. Various examples illustrate how this new approach can be used 
to find outage probability in both Poisson and non-Poisson wireless networks. 

X 

c5 ■ I. Introduction 

The spatial distribution of transmitters is critical in determining the mutual interference and 
hence the performance of a wireless network. A common and ubiquitous model for the node 
locations is the Poisson point process, where the wireless node locations are independent of 
each other. While this model offers analytical tractability, it is insufficient in many cases as it 
precludes intelligent scheduling, which leads to correlations between node locations. Another 
classical model for the node locations is the grid, but it is inflexible and too regular to model 
most wireless networks. To optimize the spatial packing of transmissions in a wireless network, 
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and obtain the maximum throughput, it is essential that smart MAC protocols are designed and 
used. In order to do so, new mathematical tools are required to analyze and understand the 
performance of general spatial networks. 

While interference and related metrics have been studied in a few non-Poisson spatial networks, 
there is no systematic approach or a general technique to analyze these problems. In this paper, 
we provide a new technique to obtain an expansion for wide range of functions of interference for 
any spatial distribution of nodes. As in a Taylor series, a better approximation of the interference 
functional can be obtained by increasing the number of terms in the series, with each additional 
term being increasingly complicated. As we shall see, this new technique can be used to analyze 
functionals of interference for a wide spectrum of MAC protocols, spatial distribution of nodes, 
and radio design choices. While in this paper we focus on treating interference at the receiver 
as noise, the techniques in this paper can be easily extended to other sophisticated interference 
mitigation techniques. 

Similar to the moments of a scalar valued random variable, the product densities [1] specify a 
stationary spatial point process. For a Poisson point process, the n-th order product densities can 
be easily obtained. This is not the case for other point processes. In practice, it is only feasible 
to obtain a few product densities (with certain reliability) by analyzing the spatial data provided. 
Therefore a non-PPP spatial distribution is usually partially specified by a few product density 
measures. The present paper uses factorial moment expansion techniques [2] to approximate an 
arbitrary function of interference using a limited number of product densities. The technique 
provided in this paper extends the rich set of results for interference and outage characterization 
available for the PPP 

A. Background and Related Work 

In a spatial wireless network, the success (or outage) probability of a typical link is an important 
metric of performance, and considerable prior work focused on obtaining the success probability 
for different spatial network models. The success probability of a typical link is equal to 



where S is the signal power of the typical transmitter, I is the interference, and N is the noise 
power. The reduced Palm measure P !o is used since we are interested in the probability of a 
"typical link", which essentially corresponds to the conditional probability for point processes. 




(1) 
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One can interpret the success probability (or other quantities) of a typical link as the success 
probability averaged over all the spatial links [3]. 

Poisson Networks: The two main reasons for the analytical tractability of a Poisson point 
process are its easy Palm characterization due to Slyvniak's theorem [1], and the knowledge of 
its probability generating functional (PGFL). 

Using the PGFL of the PPP, the Laplace transform of the interference in a PPP network is 
obtained in [4]-[8], and when the path loss model is given by ||x|| _a , a > 2, the interference 
distribution is stable with parameter 2 /a. Obtaining the outage probability in closed form using 
the distribution of interference is not always possible, but bounds on the CDF of the interference 
distribution can be obtained [9]— [1 1] which lead to bounds on the outage probability. Using the 
Laplace transform, the outage probability can be obtained in closed form for an exponential 
family of fading distributions including Rayleigh fading [5], Nakagami-m fading [12], and x 2 
and other related distributions in [13]— [16]. In [17], a novel Fourier transform-based technique 
is used to analyze the outage probability in a PPP network with arbitrary fading, but it requires 
evaluation of complex valued integrals. Using ALOHA to schedule nodes in a PPP network 
also results in a PPP transmitter set since the independence between transmitters is preserved. 
Since other MAC protocols typically induce correlation, all the above techniques implicitly or 
explicitly assume an ALOHA MAC protocol. 

Non-Poisson networks: In a PPP network, two transmitting nodes can be arbitrarily close, 
something which typically does not happen in real networks because of physical constraints 
and MAC scheduling. For example, the popular CSMA MAC protocol prevents two nearby 
nodes from transmitting at the same time. So spatial models that account for repulsion between 
transmitting nodes or other post-MAC correlations are required to model intelligent scheduling 
protocols. In [18], [19], the node locations are modeled by a PPP with a modified CSMA MAC 
protocol, resulting in a transmitter set which is modeled by a modified Matern hard-core point 
process. However, the success probability is approximated by assuming a non-homogeneous 
PPP with intensity equal to the second-order product density of the Matern hard-core process. 
In [20] the CSMA protocol is modeled by excluding interferers from a guard zone around the 
typical receiver in a Poisson model and [21], [22] analyze outage probability by using a space- 
time Poisson approximation of the transmit process and the incoming traffic. Other common 
approximation of the transmit node locations in a CSMA protocol are regular lattices [23]-[26]. 
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Clustering of nodes might occur due to environmental conditions (office spaces, gathering 
spots), be intentionally induced by the MAC protocol [27], or occur in sensor networks to 
increase the lifetime of the network [28]. In [1 1], the conditional probability generating functional 
of a Poisson clustered process is obtained, from which the outage probability is derived, and in 
[29], the outage analysis was extended to clustered networks which use intra-cluster interference 
alignment. Femtocells [30]-[32] are another example of clustered wireless networks. 

We observe that while networks that can be modeled by a Poisson point process have been 
extensively analyzed, this is not the case with networks which induce spatial correlations among 
transmit nodes. They are generally analyzed by making simplified assumptions about the spatial 
model, or by mathematical approximations in the analysis. Even in the case of Poisson networks, 
the outage analysis becomes complicated when one moves away from exponential forms of 
fading. Essentially, there is no single mathematical tool or technique that is flexible enough to 
analyze these general spatial networks. 

B. Contributions and Organization of the Paper 

The main contribution of this paper is a new mathematical tool that provides an arbitrary 
close approximation to any well-behaved functions of interference in a spatial network. We use 
this tool to provide a series approximation of the success probability in any stationary wireless 
network. Compared to a recent technique that can be used to analyze outage probability of any 
spatial network in the low-interference regime [33], the present result is far more general since 
[33] requires an asymptotically small user density. 

In this paper we use a technique first developed to evaluate the derivatives of queue functionals 
[34], [35], and later extended to obtain series expansion of functionals of Poisson point processes 
[36], [37]. The series expansion termed the factorial moment expansion was generalized to one 
dimensional non-Poisson point processes in [2], and to higher dimensional point processes in 
[38]. In [38], the FME results of [2] were extended by using measurable orderings of the points 
of the process. The main tool we use in this paper is from [2] which provides a factorial moment 
expansion of functionals of point processes. The sufficient conditions for the spatial FME to exist 
[2] are difficult to verify, and an additional contribution of this paper is the simplification of 
these sufficient conditions for functionals of interference. 

The paper is organized as follows. In Section II, we introduce the system model, and define 
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point process measures. In Section III we introduce expansion kernels and the main theorem 
which deals with sufficient conditions for the FME of interference functions. In Section IV, we 
provide the FME for outage probability, and provide various examples to illustrate the use of 
FME. 

II. System Model 

The transmitters locations are modeled by a simple, stationary, and isotropic point process 
[1], [3] $ on the plane M 2 of density A. Each transmitter x G $ is associated with a mark 
ha- that is independent and identically distributed and does not depend on the location x. The 
random variable h x may represent the transmit power, or the small-scale fading between the 
transmitter and some point on the plane. The path loss model denoted by £(x) : IR 2 — > [0, oo) 
is a non-decreasing function of ||x||, and further J B , , c £(x)dx < oo, Ve > 0, where B(x,r) 
represents a ball of radius r centered around x. The interference at y E M 2 is 

I(y,$) = 5]V(x-y). (2) 

xe<J> 

Also, let 

xe0nB(o,||z||) 

be a restriction on the interferers to the ball B(o, \\z\\). In most cases, we are interested in the 
performance of a typical transmitter and its associated receiver. We condition on the event that 
a point of the process $ is located at the origin and we consider the node at the origin as a 
typical transmitter. For example, the outage probability of a typical transmitter and its receiver 
at r(o) is given by 

Ps = P! ° Go + lWo),*) 
where S is the received signal power and N is the noise power and P !o the Palm measure [1]. 
The SINR threshold 9 for successful communication depends on the required rate, the receiver 
structure, and the coding scheme used. In this paper we shall provide a series expansion of 
functions of interference, more precisely expansions of E !o [F(l(?/, $))] for well-behaved functions 
F(x). For example, 

1) The outage probability of a typical link is equal to E !o [G(6'(l(r(o), $) + N ))], where G is 
the CCDF of S. Hence in this case F(x) = G(9x + 0N O ). 
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N +I(r(o),$) 



We now introduce a few definitions concerning point processes, and we begin by defining an 
order among the points of the process. 

Definition 1 (Measurable order). For the point process $, define an ordering of the points based 
on the distance from the origin. For a simple and stationary point process this ordering is a.s. 
unique [3] and for i,!/ G $, we denote x ^ y if \\x\\ < \\y\\. 

We now provide formal definitions of moment measures and Palm probabilities that are used 
later in this paper. Let M denote the set of finite sequences on R 2 , i.e., the set of sequences 
(j) = {xi} C R 2 such that \<p D B\ < oo for all bounded Borel B C R 2 and Xi ^ Xj, i ^ j . 
Denote by B(M) the smallest cr-algebra on M that makes the maps (f> — > \<j> n B\ measurable 
for all Borel B C R 2 . For notation simplicity we define <f)(B) = \(j)C\B\, B C R 2 . A probability 
measure P on (M, B(M)) defines a point process. Product densities characterize the distribution 
of a point process, and for many point processes are easy to characterize. 

Definition 2 (Product densities). Let Bi C R 2 , 1 < i < n. The n-th order factorial moment 
measure of a point process ($,P) is defined as 



Xl,...,X n £& 

where ^ p d represents sum over pairwise distinct tuples. The n-th order product density of a 
point process is defined in terms of the factorial moment measure by the following relation. 



So Mp (Bi, . . . , B n ) counts the mean number of n-tuples in the set Bi x B 2 x . . . x B n , with 
no two components of the n-tuple being the same, and the n-th order product density is the 
Radon-Nikodym derivative of the product measure with respect to the Lebesgue measure. For a 
stationary PPP of density A, ^"■'(xi, . . . , x n ) = \ n . This follows from the independence of node 
locations. Also for any stationary point process, the n-th order product density is a function of 
yi = x 2 - xi,. . .,y n -i = x n - xi, i.e., 



p.d 



VlP(B 1 ,...,B n )=E HxtGBt,..., 




(3) 



Q 



in \x 1 ,..., 



x n ) = Q [n '(yi, ■ ■ -,y n -i)- 
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Intuitively, the (n + l)-th order product density g( n+1 \xi, . . . ,x n ) of a stationary point process 
is proportional to the probability of finding points of the process at o,x±, . . . ,x n . In this paper 
we assume that the n-th order factorial moment measures are cx-finite, product densities of all 
orders exist for the point processes in consideration. Also, we only focus on stationary point 
processes. 

Palm probabilities are the point process counterparts of conditional probabilities of a real 
valued random variable ans are defined in terms of Campbell measures. The n-th order reduced 
Campbell measure of the point process ($,P) is a measure on (E?) n x M defined as 

A e B(M), B c (K 2 ) n . 



C in \B,A) =E 



[ l($\{x u ...,x n } e A)dQ$\xi,...,x n )) 
Jb 



Please refer to [2], [38] for a detailed description of these measures, notation and their properties. 

Definition 3 (n-fold reduced Palm measure ). The n-fold reduced Palm measure is the Radon- 
Nikodym derivative of the n-th order reduced Campbell's measure with respect to the n-th order 
factorial moment measure evaluated at {x\, . . . ,x n ). More formally, the n-fold reduced Palm 
measure is given by 

P£U(^) = dg(w)( . xA ; ) (^---^n), AeB(M). (4) 

The dot "■ " in the definition represents the variable on which the Radon-Nikodym derivative is 
defined. 

Informally, the reduced n-fold Palm measure corresponds to the law of the point process given 
that it has points at x\, . . . , x n , excluding these points. For notational simplicity we shall denote 
the 1-fold Palm measure at the origin by P !o . For a stationary point process of density A, the 
first order reduced Palm measure P !o has a simple representation [1]: 



A e B(M), 



l($x\{x}GA) 

where B is a Borel set, \B\ its Lebesgue measure, and $ x corresponds to the translation $ — x. 

In the next section, we introduce the factorial moment expansion (FME) of functionals of 
interference and provide sufficient conditions for the FME of interference functions. 
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III. Factorial Moment Expansion 

The Factorial moment expansion was introduced in [2] for point processes on the line and was 
later extended to spatial point processes in [38]. FME can be considered as a Taylor series (of 
an analytic function) kind of expansion for functionals of point processes. Similar to a Taylor 
series, the average of a functional of a point processes is represented as a finite series and an 
error term that diminishes as the number of terms increases. In a Taylor series, the terms of the 
series depend on the derivatives of the function, and in FME these derivatives are replaced by 
expansion kernels which we define below. 

A. Expansion Kernels 

We first introduce some notation. 

1) Let F(x) : IR + — > [0, oo] be a real function. Hence F(I(y, $)), is a functional from the 
space of marked point processes to real numbers, more precisely defined as 

F(l(y,$))4 F (^M(x-y)). 

xG<J> 

Hence F(l(y,$)) should be interpreted as a composition of a function F(x) and the 
interference functional, rather than F being a function of I. Similarly, F(l z (y, $)) should 
be interpreted as F(l(y, $ n B(o, \\z\\))). 

2) Adding a new point to <3> corresponds to adding a tuple (x,h x ) to the point process. But 
for notational convenience, we just represent it as $ U {x}. So when a point is added to 
the process, it implicitly means the corresponding mark (fading) is also added. 

Definition 4 (Continuity). The functional F(I(y, $))) is continuous at oo if 

lim F(I,(2/,$)) = F(I(j/,$))), 

\\z\\ — ►oo 

holds true for any simple and finite point 1 set $. 

Definition 5 (Expansion Kernels). Let F(x) : K + — > [0, oo] be a real function, the first order 
expansion kernel is defined as 

F^(l(y, $)) = F(l z (y, $) + hj(z - y)) - F(l z (y, $)), 

'A simple and finite point process is a point process for which $(£>) < oo, for \B\ < oo, and no two points coincide. 



and the n-th order expansion kernel is defined by 

FS.. A (i(y,^)) = (...(F2 ) )2 ) --0S?(i(i/^))- 

As mentioned in Remark 2, Fi"|... 2„(l(y, $)) actually means > s(I(y, $)), z'.e., 

the kernel also involves the marks of the added points and not only the points. In this sense, even 
when applied to (the point measure $ with zero mass), the kernel is random. For example, 

Fg) 22 (l( y ,$))= F (l, 2 ( 2/) $)) 

- F(l Z2 (y, $) + h z J(z 2 - y)) - F(l Z2 (y, $) + h z J( Zl - y)) 
+ F(l„(j/, $) + h Zl £(z! - y) + h Z2 £(2; 2 - j/)), z 2 ^ Zl . 

This expansion kernels can be written in a compact form [2]: 

p(B) (i(2/$))= ( Ei=o(-l) w ~ i En e {(-)}F(l*.(l/.*) + Ei 6 nM(2i-y)), 

0, otherwise, 

(5) 

where j (™) j denotes the collection of all cardinality j subsets of {1, . . . , n}. Furthermore the 
null-kernel is defined as, 

F („) (Q) = | E;=o(-l)^'En 6 {(j)} F (E j6 n - 2^,...,^, 

0, otherwise, 

Observe that F^l., Zn (0) is a random variable because of the added marks. 

B. Factorial Moment Expansion of Interference Functionals 

The following theorem deals with the FME of interference functionals, and extends Theorem 
3.2 in [2] for the reduced Palm measure and interference. Note that we state the theorem for a 
function of interference, although the main result holds for any functional of a point process. 
Let E ha . s ...^ x . denote the expectation with respect to the random variables h Xl , . . . ,h Xi . 

Theorem 1 (FME). Let F(I(?/, $)) be be such that the functional is continuous at infinity, and 
such that 

[ [ [Fg....^(l(y, ^))] I IP&V...^ - - - , ar^dan . . . da;, < oo, (7) 
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for % — 1, . . . , n + 1. Then 

n „ 

E !o F(l(y,$))=F(0) + A- 1 ^ / E^,...,^ [F« ,. (())> : ' •' (,-; r,)d.r : d.r, 

i=i "^ 2i 

+ A" 1 / / E M ,..., hin+i [F^ n+1 (l(yA))]^^ 
Jmn+i) j M 

(8) 

Proof: The proof follows the lines of Theorem 3.1 in [38]. The main difference is that while 
[38] deals with the FME of EF(l(y, $)), we focus on the reduced Palm version E !o F(l(|/, $)). 
We present the proof in the Appendix. ■ 
In short, to use FME for interference functionals, it is necessary to verify two things: 

1) the functional F(l(y, $)) is continuous at oo as in Definition 4, 

2) condition (7) is valid. 

It is difficult to check condition (7) for a general point process, and hence we now provide sim- 
plified sufficient conditions that are easy to verify. We begin by providing an easily computable 
upper bound on (7) whose finiteness can be verified. The basic idea of the upper bound on the 
expansion kernels is simple, and can be easily illustrated for the case n — 1. The first order 
expansion kernel is given by 

F^(l(y, $)) = F(l z (y, $) + h z £(z - y)) - ?{l z {y, $)). 

This can also be rewritten as 

Fi 1 \l(y,$))=hJ(z-y) [ F'(l z (y,$) + Th z £(z-y))dT, 

Jo 

where F'(x) = S^M. So if the derivative of F(x) is bounded for all x, then 

\Fm(i( y ,<s>))\<P'(x)\\M(z-y). 

Hence the first order expansion kernel is controlled by the derivative of the function F(a;), and in 
a similar manner the supremum of the n-th order expansion kernel depends on the n-th derivative 
of the function F(x), and this is made more precise in the following Theorem. 

Theorem 2. Let the function F(x) : M + — > [0, oo) be a smooth function, with derivatives up to 
order n bounded. For any 1 < k < n and 1 < < 

k 

|Fit.,, ?l (l(y,0))| < \F k \T- k l[h Xpi e(x Pi -y), (9) 
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, where \W\oo denotes the L°° norm 2 . 
Proof: See Appendix B. 

■ 

The following corollary combines all the inequalities provided by Threorem 2, and averages 
the fading. 

Corollary 6. Let 

= max{2 n - k \F k \, 0<k<n}, (10) 

then 

(ii) 

where G(x) is the CCDF of the fading random variable h. 

Proof: From the definition of J 7 *, and from Theorem 2, for any 1 < pj < n, and k < n, 

k 

\^L,My^))\<^U^n £ ( x Pi -y)- 

i=i 

Since min{l, b} min{l, a} = min{l, a, b, ab}, and using the above inequality for all k and all 
product combinations of Xi, we obtain 

,,(i(/a <t>))\ < k U mm {i, n^*, - y)} ■ 

i=i 

Since K xl ,...^yA..,xMy,<f>))]\ < E h , 1 ,...,h ;Cn [|Fg...^(l(y,0))|], and from the independence 
of h Xl , 

n 

|Eh xl ,;h, n [FW.^(I(y, 0))] | < ^ J] E ^ min I 1 ' WO* " tf)} • ( 12 ) 

i=l 

The average 

/■oo 

Eh [min{l,h a .^(x i - y)}} = / P (min {1, KA x i ~ y)} > a ) da 

Jo 

= / G Gs^)) d °- 

2 ||/|| 00 = inf {a > : L (x : \f(x)\ > a) = 0}, and L is the Lebesgue measure [39]. 



where IJ-V 



d fc F(x) 
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Alternatively, using the fact E[min{a;, y}} < min{E[x], E[y]} and from (12) we obtain the other 
bound. ■ 
The continuity of the functional at infinity is simpler and is stated in the next lemma. 

Lemma 1. If F(x) : K + — > IR + is a continuous function, then the functional F( I (?/,$)) is 
continuous at infinity as in Definition 4. 

Proof: From the continuity of F we have 

lim F V hj(x - y) = F lim V hj(x - y) 

\\Z ->0O \ '—^ I \ 2 — ><X '—^ 

\ xe<i > nB(o,||z||) / \ x64>ns(o,||2||) 

and the result follows from the monotone convergence theorem [39]. 



IV. FME of Success Probability 

In this section, we obtain a series expansion of the success probability using Theorem 1 . Each 
node x E $ is associated with a receiver r(x) at a distance R in a random direction. The success 
probability of a typical source destination link is 

P s = P !o (SINR(o, r(o)) > 6), (13) 

where, 

, , > , W sd £(i?) 
SINR(o, r(o)) - 



l(r(o),$) + N 

and N is the noise power at the receiver, and W sd is the fading 3 between the source destination 
pair. For this paper we assume W sd is independent of I(r(o)), although the case of g(l(r(o)), 
where g(x) is a random function dependent on W sd , can be dealt with in a similar manner. Hence 
the success probability is 

P s = P !o (SINR(o, r(o)) > 0) 

= P- ( g=ggS > 9 
Vl(r(o),$) + N 

= E !o F(z/N + ^l(r(o),$)), (14) 

where v = 9/£(R) and F(x) is the CCDF of the random variable W sd . The following theorem 
provides a series expansion of P s . 



3 W sd may also represent the power received by from a source by its destination at a distance x such that £(x) = 1. 
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Theorem 3. Let $ be a stationary point process of transmitters of density A such that the product 
densities are finite and such that 

X ^T* I rTmin{l,E[h]£(x fc -r(o))}^ +1) (a;i,...,x J )da;i...dx i < cx), (15) 
forl<i<n + l, where J=* = max h 1 '* d " F(v ^ vx) , < k < i\. Then, 

n „ 

P s =F(i/N ) + A" 1 V / E hxi ,... )h [F« Xi (i/N )] ^ +1 )(x 1 ,...,x J )da; 1 ,...,da; l + Err(A,n), 

(16) 

| Err(A,n)| < g±f / T min {1, E[h]^(x fc - r(o))} f? (n+2) (^i, • • • ,x n+1 )dx! . . . dx n+1 . 

(71+ 1J! y K 2(n + l) ^ 

Proof: Follows from Theorem 1 and Corollary 6. The (n + 1)! in the error term results 
from integrating over the entire domain and not a cone. ■ 
The condition (15) is a sufficient condition that is easy to verify, but not a necessary one. In 
the following subsections we provide various examples to illustrate the application of Theorem 
3. 



A. Poisson point process (PPP) 

As noted in Section II, the interference in a network with PPP distribution of nodes has been 
analyzed extensively, and the success probability has been obtained for different fading models. 
In this subsection we shall compare the existing results with the approximation obtained by the 
FME. We assume that the underlying transmitter nodes $ form a PPP of density A. For a PPP, 
the product densities are [1] 

and can easily be obtained from the independence properties. As mentioned earlier, Rayleigh 
fading, i.e., h ~ exp(l), is easy to deal with analytically using the Laplace transform of the 
interference. We begin with the FME of P s in a PPP network with Rayleigh fading. 



14 

1) Rayleigh Fading: When the small-scale fading is Rayleigh distributed, h and W sd are 
exponentially distributed. Since W sd is exponentially distributed with unit mean, from (14) the 
outage probability is given by 

P s =E !o exp(-z/N -z/l(r(o),$)), 

and hence F(x) = exp(— z/N — vx). Hence to evaluate the FME in (16), it is necessary to evaluate 
the average of the n-th order expansion kernels with respect to fading E h ...^ [Fii)...,z„(0)]. For 
n — 1, 

Eh zi [F^(0)] = E h2i exp (-i/N - vh Zl £{x - r(o))) - exp (-vN ) , 
and since h 2l is exponentially distributed, 

Similarly, the n-th order expansion kernel can be easily shown to be equal to 

Eh, 1 ,...^l F zi,...^C°)J = i , . < 17 ) 

0, otherwise, 

where 

1 



A(x) 



l + v-Hix-rio))- 1 ' 

To use Theorem 3, we have to first verify the finiteness of the integral in (15). Since F(x) = 
exp(— vx— i/N ), all its derivatives are well behaved and bounded and hence T* < exp(— z/N )(2i/)* 
for all i. The condition (15) can be simplified to checking the finiteness of 



X l F* ( [ min{l,^(x)}dx 



which is finite since min{l, i(x)} is a well behaved function. Hence using Theorem 3, it follows 
that 

P s = exp(-z/N G ) V t^- ( [ A(x)dx) + Err(A, n). 

We get the i\ in the denominator since E hzi ,...,h zn [F*i?...,2»(0)] is defined only on the simplex 
z n =4, ... ,=4 %u while the integration is over the complete domain. The error term is bounded 
by 

|Err(A,7i)| < t^—^t, exp(-^N ) ( f min {l } £(x)\ dx ] = 6(A n+1 ). (18) 
(n + ly. \J R 2 J 
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The error term tends to zero as n — > because of the (n+1)! in the denominator. We also observe 
that the error term gets smaller as the noise N increases. This is because as the noise power 
increases, the interference has smaller influence on the outage probability, and a fewer terms of 
the FME suffice for a good approximation. The exact outage analysis in a Poisson network is 
obtained in [5] using the probability generating functional, and the success probability is 

P s = exp(-z/N ) exp ( -A [ A(x)dx j = exp(-z/N ) V" ^— p- ( / A(a;)da: 
' \ ii 2 / i^o \Jr 2 

and we observe that FME provides the exact expansion. Analytically for a PPP, the difficult 

part of the FME is to obtain the expansion kernels since the number of terms in the n-th order 

expansion kernel grows exponentially with n. In Fig. 1, the FME is plotted for different orders, 

i.e., n for a = 4, N = 0, and v — 1. 

2) Nakagami-m fading: In the case of Nakagami-m fading, the CCDF of the fading is 

= r(mmx) 
v ' T(m) 

The expansion kernels do not have a simple form as in the case of Rayleigh fading. We now obtain 
the expansion kernels up to order 2 for the case of the singular path loss model £(x) = 



16 



m = 2 


27T 4 (a-4)(a + 2) i! csc ii (^) 2 -A 

? A v a 


m = 3 


2Tv' l (a-4)(a-2)(a + l) z { a + 2f csc^(^) 2 _4 
fi A f a 


m = 4 


2jr 4 (a-4)(a-2)(a + l)^(a + 2)^(3a-4)(3a + 2) ;,! csc^^) 2 4 
27c« 



TABLE I 

The third term in the FME, i.e., ^ J K2 J R2 E hzihz2 [Fii' Z2 (0)]dzidz2 for different™ 



and no noise, i.e., N = 0. For the other cases, numerical methods have to be employed to obtain 
the expansion kernels. 

V [F«(0)] = E h3i F (uh Zl £(x - r(o))) - 1. (19) 

So the first order term of the FME series is given by 

A / [E hzi F(uh Zl \\(x-r(o)\D-l]dx. 

Using simple substitutions it is equal to 

Eh 2 /" A / F {u\\x\\- a ) - 1 dx = i -^—^ aZ. 

So the FME with two terms is 

P B = l-Ai/~ * -^-^ ^ + ErrA,l. (20) 

1 [rrif 

In this case, the error is bounded by 

|Err(A,l)| < 2(z/A) 2 ( / min {1, £(x)} dx 

This error bound can be improved by using the alternative bound in Corollary 6. More terms 
in the series can be computed to have a better approximation of the error probability. For 

(2) 

example the term corresponding to E hzihz2 [Fz 1 , 22 (0)] for different m are given in Table I and the 

/ON 

term corresponding to ^ ^ [Fi 1 , Z2|Z3 (0)j are given in Table II. In Fig. 2, the approximations 
provided by FME with different number of terms is plotted for m = 2 and m = 3. 
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m = 2 


-47r e (t*-6)(c t +2) 3 csc a (^) 3 _6 


m — 3 


-47r B (a-6)(a-3)0 + l) 3 (a + 2) 3 csc 3 (^) 3 _ 6 

— S 2 — ^ ^ Q 

3a s 


m = 4 


4jr e (a-6)(a-3)(a-2)(a + l) 3 (a + 2) 3 (3a + 2) 3 csc 3 (^) 3 6 
81a." 



TABLE II 

The fourth term in the FME, i.e., ^ / R2 J R2 J R2 E hiihz2hz3 [Fi 3 1 ) , Z2 , Z3 (0)]d2idz 2 d2 3 for different m 



9 =2, a=5, m=2 =2, a=5, m=3, R 5 6=1 




Fig. 2. Illustration of the FME approximation with increasing number of terms, for a = 5 and v = 1. The left figure 
corresponds to m = 2 while the right figure corresponds to m = 3. 

5) Log-Normal Shadowing: Log-Normal shadowing is commonly used to model large scale 
fluctuations in the channel, and the CCDF of the fading is given by 

F(x) = i-lerff l0g( ^ y 

Neglecting noise, we obtain the first term of the FME as 

P s = 1 - 2Tre 4 ° 2 ^ a % 2/a \ + Err(A, 2). 

The other terms have to be obtained using numerical methods. 



B. Outage in the Matern hard core process ( CSMA ) 

The spatial distribution of the transmitters that concurrently transmit in a CSMA network 
is difficult to determine, but the transmitting set can be closely approximated by a modified 
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Matern hard-core processes [19]. We start with a Poisson point process \P of unit density. To 
each node x£$, we associate a mark m x , a uniform random variable in [0,1]. The contention 
neighborhood of a node x is the set of nodes which result in an interference power of at least 
P at x, i.e., 

A/-(x) = {yG*:h yx f(y-x)>P}. (21) 
A node x£f belongs to the final CSMA transmitting set if 

m x < m y , Vy G Af(x). 

The average number of nodes in the contention neighborhood of x G does not depend on the 
location x by the stationarity of \& and is equal to [19] 



POO 

Af = E[Af(x)} = 2tt / F (P^(x) -1 ) dx, 
Jo 



where F(x) is the CCDF of the fading distribution. The density of the modified Matern process 
$ is equal to 

/i-e^-Apy (22) 



N 

The next Lemma from [19], provides the higher order product densities of the CSMA transmitting 
set. 

Lemma 2. The n-th order product density of the modified Matern hard core process is 

Q ln) (yi,....,y n -i)=n\ [ i(0<h< ...<t n <l)f(t 1 ,...,t n )dt 1 ...dt n , 

J[o,i] n 

where 

/(*!,...,*„) =e*P<- V (-l)* J+1 t miQieJ I TTF(P£(T-y i )- 1 )dx 



[ JC{l,..,n} ,/R2 iGJ 



•niK 1 -^^-^)" 1 ))' 

wzY/i the convention y n = 0. 

Using these product densities, a simple sufficient condition for the FME to hold true is provided 
in the next Lemma. 
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Lemma 3. If J 7 * < oo, then the FME expansion holds true for the modified Matern hard core 
process and 

Kf T* / r 

I Err(A, n) \ < - _" e + xp(-AT)) \J R2 mhl {1 ' E[h]£ {x ~ r( ° ))} (1 " F ^ x ^ dx 

(23) 

Proof: It is easy to observe that 

n-l 

g^(y l ,....,y n _ 1 ) < JJl-F(P^ fc )), 

k=l 

and hence the result follows from (15). ■ 
Observe that as P — > 0, the density of the transmitter set decreases, since the contention 
set increases. In Lemma 3, we observe that the rate of decay of the error with respect to P 
depends on the behavior of F(x) at the origin. Using Theorem 3, and Lemma 3, the outage 
probability can be obtained for any fading distribution when the TXs locations are modeled by 
the CSMA-Matern process. 




C. Determinantal Point Processes (DPP) 

DPPP were introduced by O. Macchi and the points of a DPP exhibit soft repulsion. Hence 
these processes are particularly appealing for modeling node locations in structured spatial 
networks e.g., cellular networks. DPP have been used to analyze eigenvalues of random matrices, 
zeros of analytical functions, and Fermionic gases (in physics) and exhibit a rich mathematical 
structure. See [40], [41] for a good exposition of DPP. DPP is particularly suited to FME 
analysis since it is defined by its product densities. A point process is determinantal if its n-th 
order product density is given by 

p n (xi, ...,x n ) = det(K(xi, Xj))i< itj < n . 

K(x, y) is the kernel of the DPP and is assumed to be locally square integrable, a local trace 
class operator, hermitian and non-negative definite. Stationarity of a point process implies that 
the product measures are translation invariant, i.e., 

p (n) (xi +?/,..., x n + y) = p (n) (xi,...,x n ). (24) 

Hence for a stationary DPP the kernel should be of the form 

K(x, y) = K(x - y, 0) := K(x - y). 



2(1 



The average number of points in a set B C R 2 is equal to 

E[$(E)] = I p (1) (a;)dx, 
is 

and for a stationary DPP, we can observe that the density is equal to 

a^m.ko). (25) 

We have the following upper bound on the product densities of the DPP 
Lemma 4. For a stationary DPP, ifK(x) = K(— x), 

P V> ( Xl , ...,x n ) = pM(y u y n ^ 0) < m '^y}? k \ 
where y k = x k - x n . 

Proof: Follows from Fan's inequality [42]. ■ 
It also follows from the non-negative definitive nature of the kernel (A.4) that K(x) < K(0). 
We now provide a few examples of stationary DPP 

1) The Ginibre Ensemble It was proved in [43] by Ginibre that the eigenvalues of complex 
non-hermitian n x n matrices with unit Gaussian random variables, form a DPP with kernel 
(in the limit n — > oo) 

K(zi, z 2 ) = 7r _1 exp (—\\zi - z 2 \\ 2 /2) . 

The Ginibre ensemble is a stationary DPP (observe that K(zi, z 2 ) = K(z 1 —z 2 , 0)) of density 
7T _1 . The density of the point process can be modified (decreased) by changing the variance 
of the Gaussian random variables and is equal to (7rcr 2 ) -1 . See Figure 3 for a comparison 
between a PPP and a DPP 

2) The eigenvalues of A = A + ivB, where A, B belong to the Gaussian unitary ensemble, 
and v < 1, form a DPP with kernel, 

where r = (1 — v 2 )/(l + v 2 ) E [0, 1]. Observe that the modified Ginibre DPP is of density 

VWl-r 2 )). 

3) Fermionic Gas (Sine DPP) The probability distribution of fermion (n — > oo) locations on 
the real line is a DPP [40] with kernel 

sin(7r(xi - x 2 )) 
K(x 1 ,x 2 ) = -, r — . 

7T(Xi - X 2 ) 
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Fig. 3. The left figure illustrates a realization of the Ginibre DPP, while the right figure corresponds to a PPP with the same 
density. Observe that the points in the DPP seem more regular compared to the PPP. 



This has been extended to higher dimensions in [44], and for two-dimensions, the kernel 
is given by 

KM: 



2^r ' 

where r = ||xi — a? 2 1| and J\(x) = tt^ 1 JJT cos(xsin(0) — 9)d9 is the Bessel function of 
the first order. The density of the Sine DPP is K(0) = 1/2. 
The following theorem simplifies the necessary condition of the FME for the case of a DPP. 

Lemma 5. If the transmitters form a stationary and isotropic DPP with kernel K(x) and < 
oo, then the FME holds up to n terms and the error is bounded by 

-p* \n+2 / r / K 2 M\ \ n+l 

| Err(A,n)| < ^_ Ij min {1, E[h]£(x - r(o))} ^1 - ^-j dxj < oo. (26) 
Proof: From Theorem 3, it suffices to show that (15) holds true. From Lemma 4, we also 

have 



K(0) 



n—2 
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and hence (15) is 

X ~ ljr : [ f\mm{l,E[h]£(x k -r(o))}g^ +1 \x 1 ,...,x l )dx 1 ...dx l , 
<A- 1 J:; / n min I 1 . ®W(*k - r(o))} n h£^ k) dx! . . . dx t . 

^^=1 K (°) 

= A-^O) 2 ^^y^ 2 mm{l,E[h]£(x-r(o))}p( 2 )(x)dx^) . (27) 



Since for a stationary DPP, 



£ (2) (x) =K 2 (0) -K 2 (x), 



and A = K(0) > K(x). Hence g^(x) < K 2 (0) and hence (27) is always finite. We obtain the 
error bound by substituting q^ 2 \x) = K 2 (0) — K 2 (x) in (27). ■ 
Permanental point process are the counterparts of DPP which exhibit attraction between the 
points and are defined by their product densities. The product densities of a permanental point 
process are defined by the permanent of a kernel matrix, and similar to that of the DPP, the 
FME can be used to analyze outage in these processes. 

D. Real Data 

If the locations of the nodes in a real wireless system are known, can a semi-empirical formula 
of the outage probability be obtained? The answer is affirmative, and the outage can be obtained 
using the FME analysis. Given a single snapshot of the node locations, the product densities can 
be estimated using the techniques described in [45] and [46], and the FME can be obtained by 
numerical integration. Although there is an initial complexity of estimating the product measures, 
they can be reused multiple times to evaluate different functionals of the interference. 

V. Conclusions 

In this paper, we have introduced a new technique to evaluate Palm averages of functionals of 
interference in a spatial networks. A series representation of the interference function, termed as 
the factorial moment expansion (FME) was obtained. To obtain this series representation, only 
the product densities of the underlying transmitter locations are required, quantities that can be 
easily computed for many point process. A main contribution of the paper is providing simple 
sufficient conditions for the FME to hold true for functionals of interference, and computable 
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bounds on the truncation error. We have provided several examples to illustrate this procedure and 
provided bounds on error. This new technique is versatile and is limited only by the computational 
complexity and the knowledge about the underlying node distribution. 

Appendix A 
Proof of Theorem 1 

For simplicity of notation we neglect fading, and the case with fading can be easily dealt with 
since the fading is independent across nodes. From Hanish's lemma [47, Proposition 1], we have 

(Pfc*k-i)S(<^) = P K...,*„(#), v£\dx u dx n ) a.e. (28) 
Hence, using the Baccelli-Bremaud lemma [2, Lemma 3.3] applied to the simple point process 

(*,pgiU), 

Fg,...,, ! (l(2/,0))PSS..^( d ^) = F S,..^(°) 



M 



I I F^ m (l(y^))PH.^ 1 (d^)Mg! +1) (dx i+l ), 
Jm 2 J m r o,x 1 ,...,x i 



almost surely with respect to the measure Mp (dxi, ■ ■ ■ , dxi). Integrating with respect to the 

Fg,...«(l(y, 0))Pi5 1 1 , ) ..,, ! (d0)Mg(dx 1 , . . . dx,) 



measure M™\ we obtain 



!i J M 



= [ F« ..((DM^i.l.r,....,!,-,) 

+ [ [ [ FSU^^'^^Sl^iW^U (d^ + i)Mg(dx 1; ...d^ 
From [2, Proposition 2.5], we obtain 

M 2+D (dx i+ i)Mg(dxi, . . . , dx^ = M^ +1) (dxi, . . . , dz i+ i), 

and hence 

/ K 0))Pi5 1 1 , ) ..,, ! (d0)Mg(dx 1 , . . . dx,) 

J M 

= [ Fg ..((DM'i.l,-,....,!,-,) 

+ / / F^^ 1 (l(y > 0))P^„ i ^ 1 (d^)4^ 1) (da; 1 ,... > d a ; H . 1 ). 
Adding both sides of the above equation for % — 0, . . . , n, we obtain the required result. 
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Appendix B 
Proof of Theorem 2 

It can be easily seen that the difference function is also equal to 

F£U»(l(i/,0) = E(" 1 ) EILlblF ( ^niV^) +$> i h,i(* i - y) 

7>n V i=l 

where P n denotes the set of all binary tuples (b 1; . . . , b n ), bj G {0, 1} and has a cardinality 2 n . 
For notational convenience we denote I Xn (y, 3>) by /3 and h Xi £(xi — y) by 7;. So we have 

f£U.(i(v,0) = E (-i) E? = ibi F (^E b ^J • 

(bi,...,b n )e-p„ V i=l / 

Without loss of generality, we can assume Pi = i, 1 < i < k. We now partition the set V n 
into 2 n ~ k groups. Each partition consists of the n-binary string with fixed bits in the positions 

{k + l,...,n}. 

A(b k+1 , . . . , b„) = {(bi, . . . , b fc , b fe+ i, . . . , b re ), b* G {0, 1}, 1 < i < k}. 

For example with k = 2 and n = 4 we partition of the binary-4 tuples as A(0, 0), A(1,0), 
A(0, 1), and 1), with 

A(b 3 , b 4 ) = {(0, 0, b 3 , b 4 ), (1, 0, b 3 , b 4 ), (0, 1, b 3 , b 4 ), (1, 1, b 3 , b 4 )}. 

So we have 

F£U„(l(z/,0)) = E (-iF'^K,..,^)), (29) 

(«i,.. M «„_ fc )e7' n _* 

where 

(n— k 
z?+E^+fc+E b ^ 
i=l i=l 

From (29), taking the absolute value 

I4t.,,„(l(y,0))l< E |H(AK...,^_ fc ))|, (30) 

(i)i,...,D„_ fc )e'P n _ fc 

Define F k (x) = and 

tn—k k \ fc 

/3+e + e^ 1 - ^ 
i=l i=l J i=l 



25 



By a little algebra, it follows that 

■ 1 r l 



[ ■■■ [ g(9 u . . . , fc )d0i . . . d6 k = K(A(vx, . . . , v n _ k )), 
Jo Jo 



and hence 



\K(A( Vl ,...,v n - k ))\ < f ... [ \g(9 1 ,...,e k )\d9 1 ...dB k 

Jo Jo 

pi pi k 

< 1^1 / ... / \% d0i...d0 fc , 
Jo Jo i=1 

k 

= l^fcl Jj7t- 

i=l 

Substituting (31) in (30), 

k 

i4t.,x„(i(2/,0))i< e i^i 

(•ui,...,u„_ fc )e-P„_ fc i=l 

fc 

= 2- fc |^ fc |n^, 



t=i 



proving the theorem. 
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